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HOMOGENIZATION OF REACTION-DIFFUSION EQUATIONS IN FRACTURED 

POROUS MEDIA 

HERMANN DOUANLA AND JEAN LOUIS WOUKENG 


Abstract. The paper deals with the homogenization of reaction-diffusion equations with large reaction 
terms in a multi-scale porous medium. We assume that the fractures and pores are equidistributed and 
that the coefficients of the equations are periodic. Using the multi-scale convergence method, we derive 
a homogenization result whose limit problem is defined on a fixed domain and is of convection-diffusion- 
reaction type. 


1. Introduction 


Our aim is to investigate, by mean of mathematical homogenization techniques, the diffusion phenomenon 
in a multi-scale porous medium. The medium consists of a connected network made of pores and fractures 
which are equidistributed, and the diffusion process is modelled by a semilinear reaction-diffusion equation 
with a large reaction term. 

To be more precise, we consider a diffusion process modelled by the following boundary value problem: 
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X (0,r), (1.1) 

X (0,r), 


where T > 0 is a fixed real number representing the final time of the process, is a fractured porous domain 
in which the process occurred and whose structure follows in the lines below (see [22]). 

Let n be a bounded domain in {N > 3) locally located on one side of its Lipschitz continuous boundary 
i9fl. Let Y = (0,1)'^ be the unit cell in and put Y = YcU Ym where Ym and U are two disjoint open 
connected sets representing the local structure of the porous matrix and the cracks (fissures), respectively. 
We assume that a periodic repetition of Ym in is connected and has a Lipschitz continuous boundary. 
Next, we set Ym = Zg U Zp where Zp and Zg are two disjoint open connected sets representing the local 
structure of the solid part of the porous matrix and the pores, respectively. We assume that Zp and Zg have 
strictly positive Lebesgue measures and that Zg has a Lipschitz continuous boundary. The fractured porous 
medium is defined as follows. For £ > 0, we set 


Gm — UfegzN (/c-f U„) and Gc —R'^\Gm, 


and 


Gs = Ufegziv (fc-|-Zs) and Gp = Gm\Gg: 

and we define the pores space ftp = ft H e'^Gp (this include the pores crossing dfl ), the cracks space 
n® = n n eGc (this includes the cracks crossing dQ ) and the fractured porous medium as: 
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We assume that both and flp U are connected. 

This being so, the e-problem (HU is constrained as follows: 

A1 Uniform ellipticity. The matrix A(j/,t) = (fly (y,T))i<ij<Ar G is real, symmet¬ 

ric, positive definite, i.e, there exists A > 0 such that 

ll®iiII loo(rn+i) ^ a, 1 < i,j < N, 

aij{y,TKiCj > A-i ICI^ for all {y,T) G R^+\C e 

A2 Lipschitz continuity. The function g : x R x R —>■ K satisfies the following hypotheses. There 

exists C > 0 such that for any (y, r) G and u G R, 

\dugiy,T,u)\ < C 

\du 9 {y,T,ui) - dug{y,T,U 2 )\ <C\ui- M 2 I (1 + |mi| + |m 2 |)“^- 

A3 g{y, r, 0) = 0 for any (y, r) G R'^+A 
A4 Periodicity. We assume that: 

(i) g{-, •, u) G Cper(U X T) (T = (0,1)) for any u G R with Jy diVj U u) dy = 0 for all (r, u) G T x R; 

(ii) the functions lie in Lp^j.{Y x T) for all 1 < i,j < N; 

(iii) the density function p belongs to CperiY) and satisfies A“^ < p{y) < A for almost all y G R'^. 

Remark 1.1. Asa direct consequence of the periodicity and the zero mean value hypothesis for the function 
g (see precisely the first item of the hypothesis A4 above), there exists a unique i?(-, •, u) G CperiY x T) such 
that AyR{-,-,u) = g{-,-,u) and Jy R{y,T,u) dy = 0 for all r, u G R. Moreover R{-,-,u) is at least twice 
differentiable with respect to y. Furthermore, on letting G — VyR it follows from A2 and A3 that 

\G{y,T,u)\ < C\u\, \duG{y,T,u)\ < C, (1.2) 

and 

\duG{y,T,ui) - duG{y,T,u2)\ < C \ui - U 2 I (1 + |mi| + |m2|)“A (1-3) 


The motivation for the problem (HU arises from its applicability in the area of modeling of flow and 
transport in fractured porous media related to environmental and energy problems. In order to overcome 
difficulties encountered in numerical simulations in multi-scale porous media, we need to upscale such models, 
that is to find equivalent models by letting e ^ 0. This leads to model problems posed on a fixed domain 
Q = X (0,T) with suitable boundary conditions, hence relatively easy to handle numerically. 

Problem HU can also be viewed as modeling the flow of a single phase compressible fluid in a fractured 
porous medium that obey nonlinear Darcy law. In that case, is the density of the fluid, p{y) is the 
porosity of the medium while A(y, r) is the permeability of the medium. As the scale (size) of the fractures 
and that of the pores are separated (ratio of order e), we shall apply the multi-scale (or reiterated two-scale) 
convergence techniques and the framework introduced in m for the upscaling of the same problem in a 
fixed domain. 

The homogenization of parabolic equations has been widely investigated in the literature. We quote some 
works similar to ours. In m the homogenization of parabolic monotone operator in periodically perforated 
domain is considered. The problem they consider is degenerate and they use the two-scale convergence 
method to achieve their goal. In m. the authors study the homogenization of a family of parabolic equations 



div a{us,Vue) = / 


posed on a single periodically perforated domain with Dirichlet boundary conditions. In [5] is considered 
the upscaling of a convection-diffusion equation in a perforated domain made of holes periodically distributed. 
The homogenization limit for the diffusion equation with nonlinear flux condition on the boundary of a 
periodically perforated domain is studied in m- In [7] the homogenization of a semilinear parabolic equation 
in a periodically perforated domain is considered. In m the authors describe some diffusion models for 
fractured media. We also mention m where the author used the T-convergence method associated with 
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multi-scale convergence notions to get a limit law of an incompressible viscous flow in a porous medium with 
double porosity. 

Taking into account the preceding review which is very far from being exhaustive, we observe that the 
study of a problem like (EB) is relevant due to the geometry of the domain and we make use of the multi-scale 
convergence concept as the homogenization method. This leads to the following main result of the paper, 
where in the passage to the limit (as e —0), the large reaction term in the e-problem generates a convection 
term and we get a limit problem of convection-diffusion-reaction type. 


Theorem 1.1. Assume that the hypotheses AI-A4 are in place and let Ug (e > 0) be the unique solution 
to EB- Then as £ ^ 0 we have 

Ug —>• uo in LF‘{TLt), 

where uq G L'^{0,T; HqITI)) is the unique solution to 


\Zs\ p{y)dy^ ^ = div 

uq = 0 on did X (0, T) 
uq{x,Q) = u^{x) in il. 


(i(x, t)Vuo^ + divLi{x, t, Uo) — L 2 {x, t, uq) ■ Vuq — Lo{x, t, uq) in Ut 


The coefficients and operators in the theorem above are defined in Section U) 

The paper is organized as follows. The a priori estimates and compactness results are formulated and 
proved in Section[2l In Section[3l we recall the concept of multi-scale convergence and prove some preliminary 
results. Finally, Section 0] deals with the passage to the limit and the derivation of the macroscopic model 
for problem EB- 


2 . A PRIORI ESTIMATES AND COMPACTNESS RESULT 


Throughout, C denotes a generic constant independent of e that can change from one line to the next, the 
centered dot stands for the Euclidean scalar product in while the absolute value or modulus is denoted 

by I ■ I- 

With the connectedness of in mind, the space 

14 = {u e M = 0 on (2.1) 

is Hilbertian when endowed with the gradient norm, 

IkllUe = l|Vn||L2(ne) {u G K). (2.2) 

Therefore, the Lipschitzity of the function g{y,T, •) and the positivity assumption on the density function p 
readily imply (see e.g., Hire]) the existence of a unique solution Ug G LA{Q,T\Vg) fl C(0, T; L^(fl®)) to the 
problem EB- Moreover the following uniform estimates hold. 

Lemma 2.1. Assume that the hypotheses AI-A4 are satisfied. Then the following estimates hold true: 


sup llu 

i2(ns) < C'l 

(2.3) 

0<t<T 


pT 



\\'^Ueit)\\h(ne)dt <C, 

0 

(2.4) 

dug 


(2.5) 


<c, 

L2(0,T;V;') 



where C is a positive constant which does not depend on e. 

Proof. Let t G (0,T]. Multiplying the first equation in (11.11) by Ug and integrating over 11® x (0,t) yields: 


ip^^ug 


it) 


L2(n=) 


ipl 






+ 2 


oJn'^ 


A‘^\Vug{s)\^dxds = 2 


—g^ {ug{s))ug{s)dxds. (2.6) 
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But Remark 1 1.1 1 readily implies 


-9 = divC? f-, - drG (-, ■^,Ue] ■ 


£ 


£ 


£ 


which combined with (|2.6p leads to 




^ /* 2 


r2ro^+^// ^1VWe(s)| dxds < = W 

JoJn^ 

t 


L2(ns) 


f f G^{ue) ■ Vuedxds 
^ [ [ {^rG^(Ue) ■ VUe)Uedxds, 

JoJn^ 


where G‘^{ue) = G {^,-p,Ue) and drG‘^{ue) = Ue). Making use of (11.21) . the ellipticity of the 

matrix A and the boundedness of the function p, we have 

rt I- rt 


A ^ ||ME(t)||^2(oe) + 2A ^ [ [ |VwE(s)pda;ds < A ||w°||^2_, + 4G / / \ue\\'^Ue\dxds. (2.7) 

JoJn'^ ' ' JoJn^ 

For any real number <5 > 0, we have by Young’s inequality, 

4G f f ImeIIVmeI dxds < 4G^ f f \uef dxds + — f f \Vue\^ dxds. 

JoJq.<^ JoJn^ “ JoJn^ 

Choosing 5 > 0 such that j tbe inequality (IQ) yields: 

pt nt 

A”M|weWlli 2 (oe) + A-iy^ ||VwE(s)|li 2 (ne)da;ds < A||w°||^ 2 ( 0 )+4G5 dxds, 

which by means of the Gronwall’s inequality first leads to (12.3L then to (12.4F 
As for (12.51) . it follows from dm that 


.E due 

L 2 ( 0 ,T;y/) -^0 

On the one hand, (1^ and the boundedness of the matrix A imply 

rT 
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-5®(Me) 

'o 
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dt. 


[ IldivA'^VuEllyjdt < G. 

Jo 


( 2 . 8 ) 


(2.9) 


On the other hand, we have 


-5®(w) 

= sup 

/ G®(we) ■ V^pda; + / {drG^{ue) ■ VUe)ip dx 

£ 

Vy 

Jn^ Jn^ 


which by means of the Poincare’s inequality and dm yields 


-g’^iue) 


< C sup (||uE(t)||L 2 (Qe) + ||VME(t)||L 2 (a=)||V(p||L 2 (ne)) 
V' ‘/2eY,ll¥3|IVe=l 


< G (||UE(t)||L2(Qe) + ||VME(t)||L2(ne)) 
Therefore, making use of (I2.3I) - (I2.4I) and the Holder’s inequality we get 


te (o,r). 


f 

-/(w) 

Jo 

£ 


dt < G. 


VJ 


We combine (12.8p . (12.91) and (12.101) to get (12.5p . 


( 2 . 10 ) 

□ 


The next result relies on the following classical extension property (see e.g., my 
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Proposition 2.1. For any e > 0, there exists a bounded linear operator from into HQ{n) such that 
for any u £ Vg we have: 

PeU = u in (2-11) 

ll^ew||,^i(n) < C||u||u„ (2.12) 

where C is a positive constant independent of e. 

For a function u € L^{0,T; 14) we define its extension P^u as follows 

{P,u)it) = Pe{u{t)) a.e. t e (0,T), (2.13) 

and PeU € 

Bearing this in mind and owing to Proposition l2.ll (see precisely (12.121) 1. we have the following corollary. 
Corollary 2.1. Under the hypotheses of Lemma 12.11 we have the following uniform estimate 

^ ^ (2-14) 

where C > Q is a positive constant independent of e and where P^ is the extension operator defined in 
Proposition 12.11 


The next estimate requires some preliminaries. We define : ^fo(fl) —>• 14 by ReU = u|ne for u S ilo(fl) 
(where denotes the restriction of u to He). Then, R^ is continuous since 

\\Reu\\ve < l|M||i/i(n) for u G ilo(n). 

We recall that the adjoint R* : Vf ^ of R^ satisfies, for all v G Vf and ip G 7lQ(n), 

{R*v,(p) = {v,Re(p), 

where the brackets on the left hand side denote the duality pairing between the spaces and llo(12) 

while those on the right hand side denote the duality pairing between Vf and 14- It is straightforward that 

R;u = xne u for u G x (0, T)). (2.15) 


Indeed, for any ip G L'^{0,T; Hq^Q)), we have 

(i?*u, (/?) = {u,Re(p)= / / u{(p\n^)dxdt 

Jo Jn^ 

= / / Xn^iu(p)dxdt = / / {xn^ u) (fdxdt. 

Jo Jn Jo Jn 

By the way, it is worth noticing that combining (j2.15l) and Proposition 12.11 (see precisely (12.111) therein), we 
have 

R*u = xn4Peu) for u & (2.16) 

Likewise, one can easily check that, for any u G L‘^{0,T\ 14) with % G Vf), we have 


RZ 


dUe 

dt 


d{Rtu,) 

dt ■ 


We are now in a position to formulate another estimate. 

Lemma 2.2. There exists a constant C independent of e such that 

,d{PeUe) 


(p^XnO- 


dt 


< C. 


L2(o,T;ff-i(n)) 


(2.17) 


(2.18) 
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Proof. We first prove that there exists a constant C independent of e such that 

.due 


K P‘ 


dt 


< C. 


L2(0,T;ff-l(n)) 

To do this, let ip be arbitrarily fixed in L^(0,T; i7o(f2)). We have 




< 


p-^.R,v 


,dUe 


>0 


P ® 

VPVe 


dt 


at 

RsP\\l^{0.T-V^) 


L2(0,T,y/) 

< C\\Rep\\L^( 0 ,T-,v,) (see ([23) 

^ C'll‘PllL2(0,T;_Hi(n))- 
Having done this, it remains to prove that 

, f duA d{P,u,) 

But with (12.161) and (I2.17P in mind, it is easy to see that 

r,* f edUe\ ^ f dUe\ e^iXn^PsUe) ^ d{PeUe) 

= - dt -= 

and the proof is completed. 


(2.19) 


( 2 . 20 ) 


□ 


The following compactness result will be the starting point of our homogenization process. 

Theorem 2.1. Assume that the sequence (p'^Xne)e>o weakly *- converges in L°°(Xl), as e — ^ 0, to some 
real function that is different from zero almost everywhere in 11. Then the sequence {PeUe)e>o is relatively 
compact in L^(0,T; LA{TI)). 

Proof. This is a direct consequence of the convergence hypothesis on the sequence (p®xnE)e>o, Corollary [2T] 
and Lemma [221 by using O Theorem 2.3 and Remark 2.5]. □ 


3. Multi-scale convergence and preliminary convergence results 

We recall the definition and some compactness results of the multi-scale convergence theory miioiiiT]. 
We also introduce our functional setting and adapt some results of the multi-scale convergence method to 
our framework. We finally prove some preliminary convergence results needed in the homogenization process 
of the problem under consideration. We introduce the following notations: Hr = H x (0,T) and T = (0,1). 


3.1. Multi-scale convergence method. 

Definition 3.1. (i) A sequence {ue)e>o C Lf{TlT) is said to weakly multi-scale converge towards uq € 


LA{TIt X y X Z X R), and denoted 

/ UE{x,t)p{x,t,-,^,-^)dxdt- 

Iqt e e e 

for all p € L^{TtT',Cper{X x Z x T)). 


-A uo in LA{PLt), if as e ^ 0, 

uo{x, t, y, z, t)p{x, t, y, z, r) dxdtdydzdr 


Ot XY xZxT 


(3.1) 


(ii) A sequence {ue)e>Q C Lf{VlT) is said to strongly multi-scale converge towards uo € xY x Z x 

Y), and denoted uq in A (XIt), if it multi scale converges weakly to uq in A (Hr xY xZxY) 

and further satisfies 


■*^ellL2(nT) ll^ollL2(nj,xYxZxr) ^ 0. 


(3.2) 
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Remark 3.1. (i) Let u € L‘^{VtT', Cp^riYxZxT)) and define : fir —>■ M by u^{x, t) = it(x, t, f, 

for c > 0 and (x, t) G fir- Then u and u‘^ u in L'^{Qt) as e —>■ 0. We also have 

tt® —>• u in (fix)-'weak as e —>• 0, with u{x, t) = / / / u(-, •, y, z, r) dydzdr. 

J J JyxZxT 

(ii) If a sequence (ue)E>o C L‘^{VIt) multi-scale converges weakly in to some uq G xY x 

Z xT), in the sense of Definition 13.11 then (13.11) still holds for ip G C{^It\ L^^(Y x Z xT)). 

(iii) Let u G C{^t', L^^(Y x Z x T)) and define u® like in (i) above. Then u® u in L’^ip.x) as 

£ —>■ 0 . 

The following two compactness results are the cornerstones of the multi-scale convergence theory. 

Theorem 3.1. Any bounded sequence in admits a weakly multi-scale convergent subsequence. 

Let E be an ordinary sequence of real number converging to zero with e. 

Theorem 3.2. Let (ue)eeE be a bounded sequence in L^({),T]H^lLl)). There exist a subsequence E' of E 
and a triplet {uo,ui,U 2 ) G L^(0, T; ILq (D)) x L'^{VtT\L'^{T]E[p^^(Y))) x Lf‘{VtT\Lf‘(Y x T\ E[p^^{Z))) such 
that, as E' 3 s ^ 0, 


Us 

-)• 

Uo 

in 

L^{Q,T-Hl{Vt))-weak 

(3.3) 

dUe 

dxi 

w — ms 
- 

duo 

dxi 

i9ui 

+ in L^{^t) {i.<j<N). 

(3.4) 


We need to tailor Theorem 13.21 according to our needs. The functions ui and U 2 in Theorem are 
unique up to additive function of variables x, tr and x, t, y, r, respectively. It is crucial to fix the choice of 
Ml. We introduce the space 

HliYm) = {uG H^eriY) : p{y)u{y) dy = 0}, 

which is a closed subspace of Rper (^) since it is the kernel of the bounded linear functional u >—>■ fy p(y)u(y) dy 
defined on Hp^^{Y). The version of Theorem 13.21 that will be used in the sequel formulates as follows. 

Theorem 3.3. Let (ue)eeE be a bounded sequence in Lf‘({),T-,H^lLl)). There exist a subsequence E' of E 
and a triplet {uo,ui,U 2 ) G L^(0, T; ILq (D)) x L‘^{VtT', L‘^{T', HpiYm))) x L^{LIt; L^{Y x T',Hp^^{Z))) such 
that, as E' 3 s ^ 0, 


Ue 

-)• 

Uo 

in 

L\Q,T- 

Hq {LI))-weak 

(3.5) 

dUe 

dXi 

w — ms 
-)• 

duo 

dxi 

dui 

du2 

dz^ 

in L'^{LIt) {i.<j<N). 

(3.6) 


Proof. The proof is similar to the proof of [SJ Theorem 2.5]. □ 

The following weak-strong convergence result (see [THl Theorem 6] for its proof) and its corollary are 
worth recalling since they will be used in the sequel. 

Theorem 3.4. Let {ue)e^E C LP‘{LIt) and {ve)eeE C L^{Qt) be two sequences such that mq nnd 

Ve Uq in L^{Qt) with uo,vo G L^{VLt xY x Z x T). Then UeVg ™ uqVo in L^(LLt). 

Corollary 3.1. Let (ue)e^E C L'^ILIt) and {ve)e^E C L^(f2T)nL°“(D7’) be two sequences such thatu^ 

Uq and Ve - vq in L‘^{LIt) with uq, vq G L^(flT xY x Z x T). Assume further that (ve)e^E is bounded 
in L°°(LIt)- Then UeVe ™ uqVq in L'^{LIt). 
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3.2. Preliminary convergence results. We start this subsection by studying the limiting behavior of the 
sequence (xn=)e>o as e —>■ 0. To do this, we first express the characteristic function of in 17, in terms 
of those of Ym and Zg. Denoting by x% and Xp ttie characteristic functions of 17® and 17®, respectively, it 
appears that 


xl{x) = XGa (-) = XVe (-) (by y-periodicity) 

- ('-«■(?))-.(a 


Hence 


XG' (x) = 


1 - {xUx)+xl{x)) 
1 - 


(by Y and Z-periodicity). 
(a: € 17) 


(l-xx.(f))(l-xx.(j)) 

-4?)-.(?)■ 


But, XYm(‘) ® XZsi') G ^ ^hat according to {Hi) of Remark IXTl 

^ X^n Xz, m L (iZr), 

where the tensor product xy„ ® Xz. is defined by (xv„ ® Xz.) (j/> z) = xy^ {y)xz, (z), ((y, z) € 

We have proved the following result. 

Proposition 3.1. As e —0, the characteristic function xo.'^ o/17® multi-scale converges weakly in L'^IXIt) 
to XY^ ® XZp • 

We now proceed and recall properties of some functional spaces we will use. The topological dual of 
Hp(Ym) is denoted in the sequel by {H^iY^))' while L'^piY^) stands for the space of functions u € Lp^j.{Y) 
satisfying Jy p{y)u{y) dy = 0. We first recall that, since the space H^iYm) is densely embedded in L‘^(Ym), 
the following continuous embeddings hold: 

Hp{Ym) C LpiYm) C {HpiYn)) ■ 

We also recall that the topological dual of Lf{T\H^{Ym)) is (HplYm))'). This readily follows from 

the reflexiveness of the space H^iYm)- We denote the duality pairing between H^iYm) and {HpfYm))' by 
(•, •), and that of Lf{T]H^{Ym)) and Lf{T] {Hp{Ym))') by [•, •]. Thus, we have 


{u,v)=j u{y)v{y)dy 


for u e Ll{Ym) and v € HUYm), and 


[u,v]= / {u{T),v{T))dT = / / u{y,T)v{y,T)dydT 

Jo Jo Jy 


for u G L'^{T; Hp{Ym)) and v G L'^{T; {H^iYm))')- Furthermore, let VpfYm) stands for the space of functions 
u G Vpf.riY) with Jy p{y)u{y)dy = 0. Owing to the fact that the space 'Dp{Ym) is dense in Hp{Ym) the 
following result holds (see e.g. [191 Lemma 2 and Lemma 3]). 

Theorem 3.5. Let u G Vp^^iY x T) and assume that u is continuous on VpfYm) <S>L>per{T) endowed with 
the Hp(Ym)) -norm. Then u G L‘i^j.{T',{H^{Ym))'), and further 
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for all ip S VpiXm) where (•,•) denotes the duality pairing between x T) and Vpf,r(Y x 

T), whereas the right-hand side is the product of u and ip in the duality between Lpg^(T; (i?p (!"„))') and 

We now define an operator 

R-.LI^AY) ^ Ll^riY) 

U I-A XYmPU- ( 3 - 7 ) 

It is clear that i? is a non-negative and linear bounded self-adjoint operator. Using the positivity of the 
weight p we prove that the kernel of R is defined by: 

Ker{R) = {u £ L^^^fY) : u = 0 a.e. in Y^}. 

We denote by Ker{R)-^ the orthogonal of the kernel of R in Lp^^fY) while LllYm) stands for the completion 
of Ker{R)^ with respect to the norm 

\\u\\+ = WxYmP^uWLi^^i^Y) {u £ Ker{R)^). 

We denote by P the orthogonal projection from Lpg^(F) onto LQ(ym)- We recall that for u £ Lp^j.(T; Lp^^fY)) 
we define Ru and Pu by 

{Ru){t) = R{u{t)) and {Pu){t) = P{u{t)) for a.e. t G ( 0 , 1 ). 

Considered as an unbounded operator on H^lYm)), the domain of R' = XYmP^ is 

W = |m G Ll,{T-,HliYm)) : e (i7p'(Un))') j • 


We endow W with its natural norm 

I|W||W = lh||L2 (-7-.ffl(y^) -F 


du 

XY^P^ 


(m G W), 


L^AT-XHUYrn))') 


and recall an important result (see e.g., [nms] ) we will use in the sequel. 

Proposition 3.2. The operator P maps continuously W mto C([0,1]; P§(Um))> *-e-> there exists a constant 
c > 0 such that 


\\Pu\\c{[o,i]-,LAY,r.)) = sup xy^P-^u{t) 

0<T<1 


LAY) 


< c||m||vv for all u £ W. 


Moreover, 


du 


dv 

Xy,^P^,v 

— 

Xy„,P^,u 


for all u,v £ W. (3-8) 

We will also make use of the following convergence results in the forthcoming homogenization process. 

Proposition 3.3. Let p £ C^{Qt)'^P per{T) ^T>p(Ym). Let E', {ue)e&E and (uo,wi) G P^(0,P; i7g(U)) x 
Lf{VlT\L‘^{p-,Hp(Ym))) be as in Theorem AS.d[ Then 


1 /* X X t f 

lim - Ue{x,t)p{-)p{x,t,-,^)dxdt =\Zs\ [xy^ PUiix,t),p{x,t)]dx dt, 

E^s-^o e Jqe^ e e Jq^ 

where Q <\Zs\ <1 denotes the Lebesgue measure of the set Zg. 

Proof. We have 

1 /* X X t 1 /* ^ X t 

- Ue{x,f)p(-)p{x,t,-,^)dxdt = - Us{x,t)xn^{x)p(-)p{x,t,-,^)dxdt 

eJqs^ e e e Jq^ e e 

= - [ Ue{x,t)xY,r,{-)XzA^)p{-)T{x,t,-,\)dxdt 

ejn-r e £2 e £ £2 
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Bearing in mind that 


J XYrr,{y)p{y)^{x,t,y,T)dy = 0 for all (x, t,r) e x K x R), 


the same line of reasoning as in the proof of O Theorem 2.3] yields, as E B e ^ 0, 

1 

^ -/Qt 

XVmXZspuiLp dxdtdydzdr, 


r XXX X t 

■ / Ue{x,t)xY^{-)xzA-^)p{-)‘Pi.x,t,-,—)dxdt 

Jq^ e s s 

JQ, 

which concludes the proof. 


I Qt^-YxZxT 


□ 


We finally introduce the following notation-definition 

Fj = x L^^t; L^AT; HI(Y^))) x x T, H^^^AZs))), 

where Hp^^{Zs) stands for the space of functions u G Hp^^{Z) with u{z)dz = 0. We similarly define 
'Dper{Zs) and remark that the space Fj admits the following dense subspace 

= C^i^T) X (Co°°(r!T) ® VpeAT) ® Vp{Y^)) X {C^{nT) ® VpeAY X T) ® Vp,AZs)) ■ 

Moreover, Fg is a Banach space under the norm 


||(uo,Mi,'W2)||fi - l|wo||L2(o_T;ffi(f2)) + \AAL'^{O.T-,L'^{T-,Hl{Ym))) + W'^^W L'^ ■,L'^ {Y xT ■,Hl„{Z))) ■ 


4. Homogenization process and main results 
Let E be an ordinary sequence of real numbers e converging to zero with e. 


4.1. Derivation of the global limit problem. Owing to Corollary 12.II and Theorem 13.31 there exist 


(mo,'«i,M 2) G Fg 


(4.1) 


•. E' oi E such that. 

as E' 

9 e 

0 , 


PeUe 


Uq 

in 

L2(0,r; 

iLg(f2))-weak, 

d{PeUe) 

w—ms 

duo 

9m 1 

du2 

In J 

dxi 


dxi 

' dyi 

' dzi 


e ^ 0 we 

have 







PXFm ® xz, in 

L°°{n) weak- 


so that Theorem 1 2.1 1 yields 

PeUe^uo in L2(0,T;L2(H)). (4.2) 

We are now in a position to formulate the first homogenization result. 

Theorem 4.1. The triple { uo , ui , U2 ) G Fj determined above by (14.1|) is a solution to the following variational 
problem: 
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p{y)dy 


1 


f dxdt - ( 

Q.T Jq 


PXY^Ui{x,t), -^(x,t) 


G{y, T, uq) ■ VV'o dxdtdydr + 

flrxYmXT ClTXYmXT 


dxdt 

g{y, T, Mo)' 0 i dxdtdydr 


My^ '’'){VxUo + VyMi + VzU2) • (yx'fpo + VyV’i + '^z'ip2) dxdtdydzdr 


(4.3) 


ClrxYmXZs xT 


1 


{drG{y, T, uq) • {VxUo + VyUi + VzU2)) ipo dxdtdydzdr 


(IxXYmXZs xT 
, for all i' 4 ’o,'ipi,'ip 2 ) & 

Proof. Let e > 0 and let {'ipo,'tpi,'tp 2 ) G t?o°- The appropriate oscillating test function for our problem is 
defined as follows: 


X t 


X X t 


i’eix, t) = Ipo^X, t) + eifi{x, t, -,—)+£ 'tp2{x, t,{x,t) G ^ X (0, T). 


£ 


(4.4) 


Multiplying all terms in the main equation of (HU) by 4’e{x,t) and integrating over x (0,r) leads to: 

f X du f X t \ f X 

/ P{-)-^'^edxdt = - l A{-,—)VUe-V'lpedxdt + -l dxdt, 

Jql ^ (di Jnt, £ £ £ Jnt, £ £ 


or equivalently to 

/ p(-)xo"(-: "^edxdt = - f xn4-^^)A{-,^)V{PeUe)-Vifedxdt (4.5) 

J^T £ £ dt £ £^ 

+ - [ XQ'^{-,^) 9 {-,-^,PeUeW dxdt. 

£ Jqt £ £ £ £ 

We now pass to the limit in (14.5p as £1' 9 e —>■ 0 . We start with the term in the left hand side. We have 

[ '^edxdt = - f p{-)xQ.^{-,^)PeUe^^dxdt, (4.6) 

Jut £ £ £^ at Jnj, £ £ £^ at 

and we recall that 

t) = t) + £^(x t - l) + i^(x t - l) + e^^(x t - - -) + ^(x t - - 

dt ^ ^ dt dt e’ £:2^ + e Qr £’ £2'' + ^ dt ^ ’ e’ £2’ £.2^+ ^2’ £2b 


It follows from (j4.2p that 

lim f p{-)xQ.^{-,^)PeUe^^dxdt 
£ £ £2 dt 


d'tjj 

p{y)XYm {y)xz, {z)uo{x, t)^^{x, t) dxdtdydz 


QtxY xZ 




By means of Proposition 13.31 we have 
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If ,x. ,xx. d'tpi X t 


9i/>i 


p{y)XY^ {y)xz, {z)ui dxdtdydzdr 


(4.8) 


QtxF xZxT 


= l^«l / 

Jn-i 


/ i\ / iN 


dxdt 


The following limits hold: 


T f ^ \T-^ dipi, xt 

f /^\ X . dip2, XX t 

}^^oL 7’ = 0. 


(4.9) 

(4.10) 

(4.11) 


After the passage to the limit in the left hand side of (I4.11|) . we used the formula Jj- ^^dr = 0. Similar 
trivial arguments work for (14.91) and (14.101) . Thus as A' 9 e —>• 0, we have 


p(^)^0£da;d< -5> P(y)^y^ j ^^tpodxdt 

-\Zs\ [ PXYmUi{x,t),^^{x,t) 
JsIt L 


dxdt. 


(4.12) 


As regards the first term in the right hand side of (14.51) . it is classical that, as A' 9 £ —)• 0, we have 

/ A{-,^)V{PeUe)-XJi)edxdt= [ Xn^{-^ ^)A{-, ^)XJ{PeUe) ■ Vlfe dxdt 

dn|, e £ JvtT £ £ £ £ 

^ //// A{y,T){VxU^ + VyUi + V^U2){V^tfQ + Vy'4)i+V^'4)2)dxdtdydzdT. {A.li) 

nTXY„xZsXT 

Concerning the second term in the right hand side of (14.5p . we first rewrite it as follows: 

\ f X t \ f X t f X t X t 

- g{-,^,PsUe)ipedxdt = - g{-,^,PeUs)i’odxdt+ g{-,^,PeUe)'lfl{x,t,-,^)dxdt 

^ dnj, £ £ £ Jn|, £ £ dn|, £ £ £ £ 

+£ / 9{-,-^,PeUe)'if2{x,t,-,^,^)dxdt. (4.14) 

dn|, £ £^ £ £^ £'' 

It is straightforward from [B Lemma 5] that 

[ g{-,\,PeUe)ifi{x,t,-,\)dxdt = [ g{-,-^,PeUe)xn'^{-,-^)ifi{x,t,-,-^)dxdt 

jQe £ £ £ £ £ £ £ £ £ £ 


jZsl JJJ g{y,T,UQ)'tpi{x,t,y,T)dxdtdydT (4.15) 
Gt xYm xT 


as Tl' 9 e —7- 0. Likewise, it holds that 

I X t X X t 

lim £ g{-,-^,PsUe)'ip 2 {x,t,-,^,^)dxdt = 0. 

E'se^o £ £2 £ £2 £2 


(4.16) 
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It then remains to deal with the first term in the right hand side of (14.141) . We have 
\ r X t f X t 

- g{-,—,PeUe)ipodxdt = - G{-,—,PeUe)-'Vxi’odxdt 

^ Jnj, s e e e 

-[ (drG{-,^,PeUe)-VxUe]lpodxdt. (4.17) 
Jn^ \ £ £ J 


It follows from [H Lemma 5 - Remark 2] that as i?' 3 e —>■ 0, 


J (^drG{^,^,PeUs) ■VxUe'^ Ipodxdt = JJJJ {drG{y, T, uq) ■ {VxUq + WyUi + VZU 2 )) ipo dxdtdydzdr. 


flT^Y^y-Zs xT 


Likewise, 


f X t 

lim / G{-,—,PeUs)-Vxi’odxdt=\Zs 
5'3e->0 7 f 2 e e e 


(4.18) 

G(i/, r, Mo) • Vxipo dxdtdydr. (4-19) 


fir X Ym xT 


Thus, as Li' 3 e —>■ 0, we have 

If ,x t ^ \, , , 

- / g{-,^,PeUe)Wedxdt 
£ JQt, £ £ 


{drGijj, T, Mo) • (VxUq + VyUi + VzM 2 )) ipo dxdtdydzdx 


nTXYmXZsXT 


IZ, 


G{y, T, Mo) • Vx'4’0 dxdtdydr 


(4.20) 


CItxY^xT 


+ \Zs\ 

X Ynx X '7~ 

which combined with (14. 1211 - (14. 1311 concludes the proof. 


g{y, T, uo)ipi {x, t, y, r) dxdtdydr 


□ 


The second term on the left hand side of (14.31) needs further investigations. In fact, for further needs, we 
would like to rewrite it using formula (13.8p of Proposition 13.21 as follows 

dui 


PXY,„ui{x,t), ^^{x,t) 


dxdt = 


' L 


PXY^-^{x,t),'tpi{x,t) 


dxdt, 


but this requires that mi G W. 


Proposition 4.1. The function ui G Lpg^{T; Hp{Ym)) defined by (14.11) and Theorem 14.11 belongs to W. 

Proof. Let V'o = ^2 = 0 and ifi = in (|4.3L where Lp G and if G T>per{T) ® VpfVm). Using the 

arbitrariness of tp, we are led to 


, dfi, . 

PXY^Ui{X,t), -^[x,t) 


dxdt = 


XYr 


, , , 

^p-^tpdydr 

or 


YxT 


= 11 9 {y,T,uo)'fidydr - 


A{y,r){VxUo + VyUi + V zU 2 )V yfi dydzdr. 


YxT 


YmXZsXT 


But g = diVyG, and the boundedness of the matrix A implies that the linear functional 


II 9 iy,T,uo)'tpdydr - 


A{y,r){VxUo + VyUi + V^U 2 )Vyip dydzdr 


YxT 


YmXZsXT 


is continuous on 72per(T) ^VpfYm) with the Lpg^(T;i7p(17n))-norm. Thus can be extended to and element 
of Lp^riT] {Hp{Ym)y). In other words, Xy„P^ G Lp^ri^T; (HpiYm))') and the proof is completed. □ 
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Therefore the global homogenized problem of dm reads 


' {uo,ui,U2) e FJ, 

p{y)dy^ ^'^0 dxdt + 


fix’ • 

G'(j/, r, Uq) ■ Vipo dxdtdydr 


dui , X , , s 

PXY,n-^[x,t),tpi[x,t) 


dxdt 

g{y, T, Uo)tpi dxdtdydr 


fir X Ym X T 


fir X Ym X T 


1 


A{y, t){VxUo + VyMi + VzU 2 ) ■ + VyV’i + zi’ 2 ) dxdtdydzdr 


(4.21) 


flTXUmXZaXT 


1 


{drG{y, r, Uq) ’ i'^x'u-o + V^Ui + VzU 2 )) V'o dxdtdydzdr 

QTXYmXZgXT 

, for all (V’0,l/’i,i/’2) e l?g°. 

The variational problem (14.211) is termed global since it contains the macroscopic homogenized problem and 
the local problem. 

4.2. The macroscopic problem. We are now in a position to derive the equation describing the macro¬ 
scopic behavior of the e-problem dD. The variational problem (14.211) is equivalent to the following system: 




1 

Ws 

1 

Ws 


G{y, r, uq) • Vi/^o dxdtdydr 

(iTXYmXT 

A(y,r)(VxUo + Vyiti -I- VzU 2 ) • x'i’o) dxdtdydzdr 


fix X X X 


(4.22) 


{drG{y,r, ug) ■ (V^ug + VyUi + Vzit 2 )) '00 dxdtdydzdr 


fi X X Yrn Z s 'X. '~J~ 

for all 00 eCo°°(flT), 


dUl N , , S 

PXYm-^(X,t):Vl{xA) 


dxdt = JJJ g{y, r, uo) 0 i dxdtdydr 

O-TXYmXT 

A{y, r)(VxUg + VyUi + Vz'it 2 ) ■ (Vyipi) dxdtdydzdr 

CItX Ym XZsXT 

for all 01 e Co°°(nr) Vper{T) (g) Vp{Y^), 


1 


(4.23) 


and 


A{y, r){VxUg -\- Vyiti -|- VzU 2 ) • (Vz02) dxdtdydzdr = 0 

1 flTXYmXZgXT 

for all 02 e (g 'DpeAT) (g Ppe.(10n) (g Ilper(^)- 


(4.24) 


We first deal with (I4.24|) . We start with a few preliminaries. We define H^{Zs) to be the space of 
functions in H^{Zs) assuming same values on the opposites faces of Z, and satisfying u{z)dz = 0. We 
remark that if u G Hp^^{Z) with u{z)dz = 0 then its restriction to Zg (which is still denoted by u in the 
sequel) belongs to H^{Zs). Vice-versa, the extension of a function u G H^{Zs) belongs to Hp^^{Z) with 
u{z)dz = 0. We have the following result whose proof is obvious and therefore omitted. 
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Proposition 4.2. Let 1 < j < N and let {y,T) G x M. be fixed. The following microscopic local problem 
admits a solution which is uniquely defined almost everywhere in Zg. 

- xKy.r)GHl{Zg)-. 

C ^ P F) 

< / A{y,T)VzX^-y^uidz =-'^Okj / ^ (4-25) 

. for all U) € H^{Zs). 

Back to (I4.24L let '02 = ® w with ip G 'D{TLt) ® 'DperiT) 0 Dp{Ym) and w G Vper{Zs). We get 


p{x, t, y, t) dxdt 


A{y,T){VxUo + VyUi + V2U2) ■ VzUjdz 


= 0 


(IrxYmXT 

which by the arbitrariness of p gives, for fixed (x, t) G Ht and fixed {y, r) G K.^ x ffi., 


/ A{y,T)V^U2 ■V^.codz = - A{y,T){Vg;Uo + VyUi)-VzU}dz. 


(4.26) 


JZs JZs 

By inspection of the microscopic problems (14.251) and (I4.26|) it appears by the uniqueness of the solution to 
(I4.25P that for almost all (x, t, y, r) fixed in Ltr x x R. 


/ \ f ^^0 / \ dU\ , 1 / \ • ry 

U 2 (a:, 0 y,r) = 2 ^ —(x,t) + —(x, 0 y,r)j x 0 ?/,r) a.e., m Zg. 

For further needs, we introduce a notation. We define the matrix V^x by 

Then we can write in short for almost all (x, t, y, r, z) G Ht x R^ x R x R^, 


(4.27) 


VzM 2 = VzX • (VxMO + VyUl). 


(4.28) 


We can now proceed and look at the mesoscopic scale. Let (x, t) G Ht and (r, ^) G R x R'^ be freely fixed 
and let 7r(x,t,r, ^) be defined by the mesoscopic cell problem: 


' 7ri(x,0r,^) G W 

PXY^^p-,fji dxdt- 

OT 




JJ g{y,r,r)'ipi dydi 


YmXZsXT 


(4.29) 


YmXT 
for all 01 G W, 


Where I stands for the identity matrix. For the sake of simplicity, we put A = A{I + VzX) dz. The 
following Proposition addresses the question of existence and uniqueness of the solution to the variational 
problem (14.291) . 
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Proposition 4.3. 


The following local variational problem admits a solution which is uniquely defined on 


' Tri{x, t,r,^) e W 

r stti 

YmXT 

g{y,T,r)'ilJidydT + JJ A{y,T)^ ■ Vyipidydr 

YmXT Y„xT 

for all 'tpi G W. 



- 11 My, T)Vy7ri • Vy-tpi dydr 


(4.30) 


Proof. It is clear from the boundedness of the bilinear form on the left hand side, and the boundedness of 
the linear form on the right hand side of (|4.30|) that (|4.30|) admits at least one solution in W. As for the 
question of uniqueness, let tti, 9i be two solutions to (14.301) . Then (^i = m — solves (14.301) with zero right 
hand side. This yields: 

- // = 0 . 

YmXT 


PXY„ 


5r ’ 


But formula (13.81) of Propositionimplies 


PXY^ 



= 0 . 


We are left with 

JJ A(i/,T)VyCi • V^Cidydr = 0, 

YmXT 

which by the uniform ellipticity of the homogenized matrix A yields Vyfi = 0 a.e. in Ym x T. Therefore 
there exists a function h depending only on r such that fiiy^r) = /i(t) for almost every {y,T) € Ym x T. 
But, since G Hp{Ym)), we have 


0 = / p{y)^i{y,T)dy = h{T) / piy)dy a.e. in T. 


'Y„ 


(4.31) 


Thus h = 0 since Jy p{y) dy ^ 0. Therefore = 0 almost every where in Ym x T. 


□ 


Taking in particular r = wo(a:, t) and ^ = Vwo(x, t) with (x, t) arbitrarily chosen in fix and then choosing 
in (14.291) the particular test functions ifi = (p{x,t)vi, with p G and Vi G i'Dperi.T) ®T)p{Ym)), and 

finally comparing the resulting equation with (14.231) . it follows by means of Proposition 14.31 (bear in mind 
that 'Dper{T) ®T)p{Ym) is dense in W), that for almost every {x,t) S fir we have 


ui{x,t) = TTi{x,t,uo{x,t),'^xUo{x,t)) (4.32) 

almost everywhere on Ym x 'T. The linearity of the problem (14.301) suggests a more flexible expression of its 
solution TTi. We formulate the following variational problems 


' uji{x, t,r) eW 


duji ' 
PXYm-^,W 


for all tp G W. 


— 11 A{y,T)VyUJi-Vyiidydr = 


g{y,T,r)ipdyd7 


(4.33) 


YmXT 


YmXT 
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and 

' 0{x,t) = 6 (W)^ 

- JJ ^(2/,T)Vy6»i ■ Vyi/idydr = 1^^ JJ Aik{y,T)^dydT (4.34) 

YmXT 

for all i/j £ W, 

and leave to the reader to check that they admits solutions that are uniquely defined on Ym x r and satisfy 

T^iix,t,r,Oiy,T) = S{x,t,y,T) ■ i + uJi{x,t,y,T,r). (4.35) 

Hence, the same lines of reasoning as above yields: 

ui{x,t,y,T) = 0{x,t,y,T) ■ Va:Uo{x,t) + u;i{x,t, V, T, uo(x, t))■ (4.36) 

We are now in a position to formulate the strong form of the macroscopic variational problem (14.221) . 
Substituting in (j4.22l) the expression of VzU 2 obtained in (j4.28p . we have: 


00, , 


p{y)dy ) I ^ “ 


G{y, T, uo) • V '00 dxdtdydi 


^T^ym'XT 


1 


~ jy] JJJ ^(y^'^)(y=>:'^o + '^yUi) ■ i^x'ipo)dxdtdydT 

0,T X Ym X 7” 


1 


SlTXYmXT 

for all V'oeCo“(HT). 


drG{y, T, Uo) • f fy (^ + X7zx)dz\ ■ (V^zUo + V^ui) 


ipQ dxdtdydr 


We put B{y, r) = fz (d + X^zX(y> y, z))dz {{y, r) £ x R) and use (14.361) to get: 


(4.37) 


'Y„ 


p{y)dy ) y ^'*^0 = “ 




drG{y,T, Uo) ■ Va,uo ) i^o dxdtdydr 
Y„,xT ) 


- 7 ^ [ AVz;Uo ■ Va;i’o dxdt --Y I 

\^s\ Jqt I^sI d KJJy^xT 

Qt 


1 

Ws 

1 

Ws 


A{y, r)VyWi dydr j ■ V^^’o dxdt 

I { (// r, Uo) ■ {B{y, t){I + Vyd)^ dydr^ ■ Va;lto|’ V'O dxdt 

six' 

I (^JJ drG{y,T,uo) ■ (^B{y,T){I+ Vy0)-VyUJi^ dydr'^ tpodxdt 


for all V'oeCo“(HT), 


where 


A{x,t) = I A{y,T){I+ Vy0{xA,y,T))dydT {{x,t)£^T)- 
'Y^xT 


(4.38) 
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Setting 


Li{x,t,r) 

L2{x,t,r) 

L'i{x,t,r) 


YmXT 


My, T)\/yU;i{x, t, y, T, r) dydr, 


dydr, 


JJ ^ \Zs\drG{y,T,r) + drG{y,T,r) ■ (^B{y,T){I+ Vy0{x,t,y,T))^ 
jj drG{y,T,r) ■ {B{y,T){I + Vy6»(a;, t, y, t)) • VyWi(a;,t,j/,T,r)) dydr, 


'YmXT 

we are led to the following result. 

Theorem 4.2. The function uq determined hy (sm and solution to the variational problem (14.221) . is the 
unique solution to the following boundary value problem: 

jZsj p{y)dy^ = div (^A{x, t)'^ uo^ + div Li{x,t,uo) — L 2 {x,t,uo) ■ Vmq - L3{x,t,uo) in 

uo = 0 on dil X (0,T) 


Uq{x, 0) = u°(j;) 


n. 


(4.39) 

Proof. The claim that uq solves the problem (14.391) has been proved above and the uniqueness of the solution 
to (I4.39|) follows from the fact that the functions Li{x,t, ■) (1 ^ ^ Lipschitz. This can be proved 

by mimicking the reasoning in [3] □ 

We can now formulate the homogenization result for problem dm. 


Theorem 4.3. Assuming that the hypotheses Al-Af are in place and letting Ue (e > 0) be the unique 
solution to dm, we have, as e — >■ 0 , 

Ug —>• Uo in LF‘{TLt), 

where uq G L'^{0,T; is the unique solution to (I4.39p . 
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